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ABSTRACT 

The design of sampled analog filters using Z-transform 
techniques developed in digital filter theory is studied. 
Both recursive and non-recursive algorithms are implemented 
and analyzed. Some deviations from theoretical frequency 
responses are noted. An analog Fourier transformer is 
constructed by cascading a recursive comb filter and a 
non-recursive bandpass filter. 
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I. 



INTRODUCTION 



Whether one is extracting information, such as separating 
signals from noise in a communications system, detecting 
targets in a radar system, or preparting signals for trans- 
mission or display, signal processing is a major portion 
of system design. This processing may be accomplished by 
several means, ranging from simple RLC circuits to complex 
digital computers. Several factors will influence the 
choice of processing scheme, among these are nature of the 
signal and noise, desired results, weight, power consum.ption , 
cost, and convenience. Each scheme has its own advantages 
and liabilities. 

One such type of signal processing scheme is the discrete 
time or sampled signal processing. This method utilizes a 
signal which is sampled at regular intervals. The sampled 
values are combined with previous and subsequent samples by 
some algorithm to produce the processed output. Implementa- 
tions of these algorithms are referred to as sampled filters. 
Basic Fourier analysis indicates that periodic sampling in 
time makes the filter response periodic in frequency [1]. 
Busiginies and Dishal report a sampled filter, developed by 
Clark, in 1949 [2] . The repetitious nature of the frequency 
response became useful with the development of radar systems. 

A pulsed radar signal has a discrete spectrum, with 
spectral lines spaced from the carrier frequency at m.ultiples 
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of the pulse repetition frequency (PRF) [3] . The signal 
to noise ratio of a received signal could be improved by 
passing the signal with the noise through a filter with a 
series of narrow passbands also spaced at the PRF. The 
signal components would all be passed by the filter while 
only the noise components close to the signal frequencies 
would pass. The remaining noise is attenuated by the filter. 
Such a filter is called a "comb" filter due to the appearance 
of the passbands as teeth of a comb. Figure 1.1 illustrates 
the use of a comb filter. 
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FIGURE 1.1. Signal to Noise Ratio Improvement 
by Use of a Comb Filter 
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While it is theoretically possible to build a comb filter 
out of many individual bandpass filters, maintaining the 
proper harmonic spacing between them is difficult. A 
sampled filter will have a frequency response which is 
periodic at a constant frequency defined as the repetition 
period, f . A low pass sampled filter, as shov;n in Figure 
1.2, will be repeated as a bandpass filter centered at 
multiples of f , hence the sampled low pass filter forms a 
comb filter. 



Continuous 
Low Pass 
Filter 




Sampled 
Low Pass 
Filter 




FIGURE 1.2. Sampled Low Pass Filter 

For a given continuous low pass filter with a cut-off 
frequency of f , the resulting derived sairtpled filter will 
have a bandpass bandwidth of 2f^. Since a narrower ’tooth' 
will reject more noise, the cut-off frequency of the sampled 
filter should be a small fraction of the repetition period. 

Another application for comb filters is found in radar 
systems which use Doppler shift to detect moving targets. 



14 



Moving Target Indicator (MTI) radar uses a repetitious 
filter to suppress the 'clutter' returns from stationary 
targets while allowing shifted returns from moving targets 
to pass. This canceller filter can be implemented by a 
sampled high pass filter, as illustrated in Figure 1.3. 
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FIGURE 1.3 Sampled High Pass Filter as a Canceller 

The sampled implementation of comb filters may be 
accomplished by a recursive technique where the previous 
inputs and outputs are weighted and fed-back to the input. 
The results of recursive algorithms are treated in great 
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detail by Gold and Rader in_£eJ.eiie iice 1. The values of 
the weighting coefficients determine the filter cut-off 
frequency. Chapter II of this study looks at the details of 
the relationships between the coefficients and the filter 
response. 

Current literature suggests areas other than radar where 
comb filters are applicable. Arndt suggests the use of a 
comb filters for TV signal enhancement in noisy channels [4] . 
Illetschko and Mavrodieva utilize comb filters in color TV 
processing of chrominance, luminance, and brightness compo- 
nents [5], [6]. Uses for comb filters in digital controllers 
are suggested by White [7]. Other areas of application will 
no doubt emerge with time . 

Since a comb filter effectively 'samples' the frequency 
spectrum of the input signal with its viniformly spaced 
passbands, information of the signal's spectral content is 
available in each comb tooth. A tunable bandpass filter used 
in conjunction with a comb filter would make it possible to 
isolate the spectral content of a single tooth. Observation 
of the behavior of that signal spectral component in 'real 
time' would be possible. Webb indicates that cases where 
short time variations in spectral components is of interest 
arise in the aerospace discipline [8] . 

As will be shown in Chapter III, a txinable bandpass 
filter can be implemented by sampled techniques. By taking 
the impulse response of the desired continuous bandpass 
filter response, sampling it at an interval equal to the 
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delay between filter outputs, and setting coefficients for 
weighting the signal samples proportional to the impulse 
response sample values, one may approximate the continuous 
filter response. Such a technique is non-recursive. Filters 
implemented by the non-recursive method are called transversal 
or finite impulse response (FIR) filters. 

The algorithms used in sampled signal processing tech- 
niques require the operations of delay (storage) , multiplica- 
tion (weighting) , and addition. The sampled signal values 
may be stored digitally or by analog means. The weighting 
coefficients must also be stored by one of these means. 

Either analog or digital methods may be used to implement 
the arithmetic operations. Until the recent development of 
charge transport devices (CTD) , only sampled digital signals 
could be stored. Since digital storage and logic hardware 
has been increasingly available and inexpensive, a great 
deal of study has been done on the digital implementation 
of sampled filters, or simply, digital filters. Now analog 
'shift registers' are available for implementing sampled 
analog filters. This study will show that the various 
Z-transform design methods developed for digital filter theory 
may be equally well applied to sampled analog signal proces- 
sors. As a demonstration, a recursive comb filter and a " 
non-recursive bandpass filter will be implemented and 
cascaded to form, an analog Fourier transformer. 
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II. SAMPLED ANALOG RECURSIVE FILTER 



A. THEORY 

The block diagram of a recursive sampled filter is 
illustrated in Figure 2.1. Since this diagram has two 
delays in that particular configuration, it is referred to 
as a second-order recursive filter. 




FIGURE 2.1 Second-Order Recursive Filter 
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Analysis of such a system, performed in Appendix A, indicates 
that the transfer function, H(z) , for this configuration is 



given by 




ag + a^z 



-1 




-2 




-1 




-2 



where z is defined as and T is the time delay 



For the special case where a2 = b2 = 0 , the system becomes 
a first-order filter. Higher order filters may be obtained 
by paralleling, cascading, or both, first and second-order 
sections . 

The frequency response of a sampled filter can be 
determined from its transfer function, H(z), by using the 
definition of z. 



Analysis performed in Appendix B indicates that the delay 
terms introduce sine and cosine variations in the frequency 
response, H(ja)), which yield a periodic nature. The 'period' 
of the repetitive response is defined as f , where 



H(jw) = H(z) 



jtoT 



z 



e 
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Two major design methods are available for determining 
the transfer function coefficients: the standard Z-transform 

and the bilinear Z-transform. Both methods will be studied. 

1 . Standard Z-Transform 

The basis for the standard Z-transform comes from 
the impulse response of the desired filter response, h(t) . 

The impulse response is sampled at the regular interval, T, 
one time delay, which corresponds to z ^ in the Z-domain. 

The Z-transform is simply a power series in z~^ where the 
coefficients of the series are the sampled impulse response 
values, or 



K(z) = r h(kT) z ^ . 
k=0 

For many impulse responses of interest this series can be 
expressed in closed form. Tables listing these closed forms 
are called the Z-transform tables. Since the designer often 
knows the Laplace transfer function, H(s), of the desired 
filter, some tables are constructed to allow one to bypass 
the actual impulse response and go directly from H(s) to 
H (z) . 

a. First-Order Low Pass Filter 

For a simple RC low pass filter, the impulse 
response is known, 

-Wet 

h(t) = w^ e , 
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where ^ cut-off frequency. This impulse response 

has a Z-transform, 



w. 



H(z) = 



1 - e 



-w T -1 * 

c z 



The 3 dB corner frequency, co^, is found in Appendix B to be 
related to the feedback coefficient, bj^, by the following 
equation. 




FIGURE 2.2 0 ) vs b, 

c 1 

For a normalized transfer function, ag = 1 and a^^ = 0. 
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b. Second-Order Low Pass Filter 



The impulse response for a second-order LC 
low pass filter is, 

h(t) = w sin 0 ) t , 

C w 

where w = (LC) ^ . The resulting Z-transform is, 
c 



H(z) = 



w sin w T z 
c c 



-1 



1 - 2cos w^T z ^ + z ^ 



For a normalized transfer function. 



= 



a, = 



a~ = 



0 , 

1 / 
0 , 



b, = 



-2cos w^T , 



bo = 1 . 



(D^ and b^ are related by the following equation, which is 
graphed in Figure 2.3, 




_ 1 _ 

2tt 



cos 



-1 
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FIGURE 2.3 vs b, 

c 1 



c. First-Order High Pass Filter 

The impulse response for an RC high pass filter 
is found in Appendix B to be, 



h (t) = 6 (t) -we 

c 



-w^t 



The resulting Z-transform is. 



H(z) 




1 




T 



-1 

z 
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By inspection, the coefficients are. 



^0 







a 



1 




/ 



b 



1 




Another first-order high pass filter can be 
constructed from the low pass filter by changing the 
polarity of b^^. 



H(z) = 



■“c'^ 
e z 



1 + 



The polarity change has the effect of shifting the response 
by ^ , so the new corner frequency is ^ - co , where w 

J. J. w w 

and are related by Figure 2.2. 

d. Second-Order High Pass Filter 

The impulse response for an LC high pass filter 

is given by. 



h(t) = (t) - sin u t , 

w 

2 1 

where co = ^ 1 =: • This yields the following second-order 

C J_iL. 

Z-transform, 
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H(Z) 



-2 



1 - 2cos 0 ) T z 
c 



+ z 



thus 




-2cos 0 ) T - w sin w T , 
c c c 



a2 - 1 / 



b, = -2cos u) T 
1 c 



2 . Bilinear Z-Transform 

For the bilinear Z-transf orm, the Laplace transfer 
function, H(s), is used to determine the Z-transform 
directly by using the following approximation for s as 
described in Reference 1, 



This formula has the effect of warping the frequency of the 
desired filter response, H(jto), so that. 






bilinear 




Appendix B shows the precise derivation of this relation. 



a. First-Order Low Pass Filter 



The transfer function for a single pole low pass 

filter is, 



H(s) = 






S + (i). 



The result of substitution for s is the bilinear Z-transform, 



(jj T 0) T 

c ^ c 



-1 



H(z) = 



o)T + 2 o)T + 2 
c c 

0) T - 2 - 

1 1TTT2 ^ 

c 



A normalized transfer function would have coefficients. 



1 , 

1 , 

0) T - 2 
c 

u) T + 2 * 
c 

The relationship between the bilinear Z corner frequency, 
0 ) , , and b, is given by the following equation, which is 
graphed in Figure 2.4, 




“c‘ 1 -1 ^l"^^ 

o)_ ir b, - 1' 
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FIGURE 2.4 w , VS b, 
c 1 



Since Sq = a^, the bilinear Z low pass filter has a zero for 

the frequency , which means the low pass response goes 

s 

all the way to zero. 

b. Second-Order Low Pass Filter 

The transfer function for an LC low pass filter 



is, 

= -rf-i ■ 

s + 

c 

—I' 

where o)_ = (LC) ^ . The resulting normalized bilinear 
Z-transform is given by 



H(s) = 



l + 2z^+z^ 



2w V - 8 _2 

1 + ( — ^ ) z -^ + z 



2 2 
0) T + 4 
c 



The corner frequency, , is dependent on b^^ by the following 

relation, graphed in Figure 2.5, 
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FIGURE 2.5 w , vs b, 
c' 1 

c. High Pass Filters 

Analysis performed in Appendix B shows that 
for the bilinear Z-transform high pass filters, the same 
form, and coefficient values with the exception of the 
polarity of a^^ as in the low pass case are used. The corner 
frequencies are found using the same design curves as for 
the low pass case. 

d. Advantages of Bilinear Z-Transform 

The bilinear Z-transform design method has 
advantages over the standard Z-transform.. Given a Laplace 
transfer function, H(s), the Z-transformed transfer fvinction, 
H(z) , can be obtained directly. In addition, the bilinear 
Z-transform places a zero to force the stop band attenuation 
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dovm . The position of this zero may be shifted to transform 
a high pass filter into a low pass filter, and vice versa, 
by reversing the polarity of the coefficient, a^. The a^ 
coefficients of a bilinear Z-transform are always integral 
multiples of one another, reducing design complexity. These 
advantages occur at a price. The frequency scaling of the 
bilinear Z-transform response is not linear with the frequency 
response of H(s) , but rather follows a tangent function. As 
a result of this frequency warping, the bilinear Z-transform 
frequency response will only approximate that of H(s) for 
frequencies where. 



The response cut-off frequency may have to be compensated, 
pre-warped, in order to occur at the desired value. 

3 . Sampling Effects 

For either of the design methods described above, 

the frequency response will repeat at multiples of the 

repetition period, f , so that the response shape for 

s 

0 < f < hi is mirrored from hi ^ f<= repeated 

from kf < f < (k+1) f for k = 1,2,3,... . Figure 2.6 
illustrates the repeated response of an arbitrary filter 
response. 
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FIGURE 2.6 Repetitious Frequency Response 

This repetitious response effect is equivalent to considering 
an input at an arbitrary frequency, f, such that 

f-kf=f , 0<f< i$f^ . 

s a — a — ^ s 

th 

That is, the input is 'heterodyned* with the particular 
multiple of the repetition period such that the difference 
frequency, f , is less than J^f . The filter frequency response 

3. S 

H(f_) is then applied to the input. Therefore, for any input 

CL 

signal frequency, only the filter response from 0 £ 
is used. If the input frequency meets the Nyquist criterion, 
i.e. is less than hf , the multiple of the repetition period 
used to heterodyne the signal, k, is zero. For frequencies 
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beyond the Nyquist bandwidth of half the repetition period, 
k is some positive integer. The ambiguity of k is an 
aliasing effect. For k = 0, there is no ambiguity, thus no 
aliasing. 

B. EXPERIMENT 

1. System Configuration 

A first-order recursive filter was implemented using 
a Reticon Corporation Serial Analog Delay (SAB-100) device 
as the delay element. The transfer function coefficients 
were obtained by potentiometer voltage dividers . The 
arithmetic operations were implemented using operational 
amplifiers with feedback. Figure 2.7 illustrates the system. 



'in 




'out 



FIGURE 2.7 First-Order Recursive Filter 
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A variable feedback resistance was installed on the Reticon 
SAD-100 driver card output amplifier so that the gain of 
the delay circuit, -a, could be adjusted. This factor, a, 
is important to the setting of the coefficient potentiometers. 

Circuit analysis of the coefficient potentiometers, 
performed in Appendix C, indicates that the maximum coeffi- 
cient values for a^^ and b^^ is a. In addition, since the 
input impedance of the summing circuit is finite, the 
coefficients are determined by a loaded voltage divider 
relationship graphed in Figure 2 . 8 for the particular loading 
used in the system. 




FIGURE 2.8 Coefficient Calibration Curve 
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The coefficient calibration curve may be experimentally 
checked by observing the filter response to a unit step 
input. Appendix D derives the expressions for the coeffi- 
cients in terms of the differences between successive filter 
outputs. Figure 2.9 shows the step response for a typical 
first-order lov; pass filter. The coefficients are calculated 
by the following relationships, 



a 



0 







for k 



3,4,5, . 



This technique is extended to second-order filters by the 
use of the following relationships. 



a, = 



di , 

^3^ - 

d3^ - 2d^d3d^ t d^^d3 t d^d/ - d^d3d3 

s' - SS 



b, = 



^2*^5 ~ ^3S 

s' - ss 



s - ss 

^3' ■ ‘^2‘^4 



Appendix E describes the use of the Tektronix 31 programable 
calculator to solve these equations. 
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For first-order standard Z-transform filters of the 



form. 



H(z) = i Y • 

1 + b^z 

another method for experimentally determining is available. 
As discussed in Appendix B, by taking the ratio between the 
response minimum and maximum, y, b^ is determined by, 

^l' 1 + y 



which is graphed in Figure 2.9. 




FIGURE 2.9 bj^ vs y 
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2 . Experimental Frequency Response 

Several types of filters were implemented to verify 
the design theory and observe the characteristics of tne 
repetitious response. For an implementation using a delay 
time, T, theory indicates that the filter response should 
repeat exactly at multiples of f^, where. 




Since the SAD-100 is a 98-bit register, the delay time is 
98 times the clock period, T^, 

T = 98T . 

c 

The repetition period would then be the clock frequency, 
f^, divided by 98. Since a new input is sampled every clock 
cycle, the rate of sampling is f . The sampling is done at 
the input gate of the SAD where a transistor is turned on 
by one of the two clock phases for half of its period. 

Since the two phase clock operates at the 'on' time, 

i.e. the sam.pling pulse width, is T . The finite pulse 

width causes an envelope over the sampled signal spectrxam. 

s in X • 

The envelope follows a shape. This envelope has its 

X 

first null at f = = f^ . Since the clock frequency is 

c 

N times the repetition period, where N is the number of 
delay bits, N frequency response cycles will occur before 
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the first sampling envelope null. Because of aliasing from 
the sampling rate, f^, only those responses for frequencies 
less than or equal to i.e. the first JjN, are useful. 

In this study, only the first twelve cycles- are observed, 
therefore the minimum value of the envelope was. 



, 12 , 

sxn tt(^) 

IT (gg) 



0.976 , 



which is sufficiently close to unity to be neglected. 

Standard Z-transform filters were implemented with. 



ag = 1 , a^ = 0 , b^ = -0.3, -0.7 low pass, 

b^ = 0.3, 0.7 high pass. 



Bilinear Z-transform filters were implemented with. 



ao=l, a^=l, b^= -0.7, 0.0, 0.3, 0.7 low pass, 



a„ = 1 , a, = -1 , b, = -0.3, -0.7 



high pass, 



Two clock frequencies were used. 



f = 400 kHz, 1 MHz . 

Figures 2.10a through 2.19a illustrate the first response 
period of the various filters implemented. The solid curve 
represents the calculated response while the points are 
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the observed response. Figures 2.10b through 2.19c illus- 
trate several response cycles of the filters. 

Some significant deviations from the theoretical 
responses were noticed, especially for the larger values of 
bj^. The responses are not mirror images of one another. 

The reason for this 'skewness' is not clear at present. One 
possible explanation may center on a frequency varying 
coefficient, such as a^. The values of Uq and a^ must stay 
equal for a zero to appear at multiples of f . If one of 
the values changes differently with frequency than the 
other, say through a frequency characteristic of the delay 
device, the zero would not be properly formed, perhaps 
shifted up or down frequency from kf^. The other major 
deviation is the change of the maximum amplitudes of the 
responses with frequency rather than staying uniform as 
theory predicts. 

For filters with small b^^, the trend is for decreasing 
amplitude v;ith increasing frequency, particulary beyond 40 kHz. 
This effect is probably due to the gain roll-off of the 
operational amplifiers in the summers. Part of this decrease 
could be due to the sampling envelope, though it is always 
greater than 0.97. 

For bilinear Z filters with b^^ = ±0.7, the maxima 
of the frequency responses increase with frequency, as 
shown in Figures 2.11 and 2.14. This trend may also be due 
to deterioration of the zero from aQ and a^ varying more 
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with increased frequency. The effects of frequency varying 
coefficients need to be pursued. 

Measurements of corner frequencies, which are 
STOimnarized in Table I, indicated good agreement with the 
design theory. 
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FIGURE 2.10a Theoretical and Observed Frequency Response 
of Standard Z Low Pass Filter; = -0.3 
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FIGURE 2.10b Theoretical and Observed Frequency Response 
of Standard Z Low Pass Filter; bj^ = -0.3 
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FIGURE 2.10c Theoretical and Observed Frequency Response 
of Standard Z Low Pass Filter; = -0.3 
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FIGURE 2.11a 



Theoretical and Observed Frequency Response 
of Standard Z Low Pass Filter; = -0.7 
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FIGURE 2.11b Theoretical and Observed Frequency Response 
of Standard Z Low Pass Filter; bj^ = -0.7 
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FIGURE 2.11c Theoretical and Observed Frequency Response 
of Standard Z Lov7 Pass Filter; = -0.7 
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FIGURE 2.12a Theoretical and Observed Frequency Response 
of Standard Z High Pass Filter; = 0.3 
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FIGURE 2.12b Theoretical and Observed Frequency Response 
of Standard Z High Pass Filter; b^^ = 0.3 
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FIGURE 2.12c Theoretical and Observed Frequency Response 
of Standard Z High Pass Filter; = 0.3 
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FIGURE 2.13a Theoretical and Observed Frequency Response 
of Standard Z High Pass Filter; = 0.7 
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FIGURE 2.13b 



Theoretical and Observed Frequency Response 
of Standard Z High Pass Filter; h-^ = 0.7 
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FIGURE 2.13c Theoretical and Observed Frequency Response 
of Standard 7. High Pass Filter; = 0.7 
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figure 2.14a Theoretical and Observed Frequency Response 
of Bilinear Z Low Pass Filter; = -0.7 
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FIGURE 2.14b 



Theoretical and Observed Frequency Response 
of Bilinear Z Low Pass Filter; b^ = -0.7 
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FIGURE 2.14c Theoretical and Observed Frequency Response 
of Bilinear Z Low Pass Filter; = -0.7 
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FIGURE 2.15a Theoretical and Observed Frequency Response 
of Bilinear Low Pass Filter? = 0.0 
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FIGURE 2.15b Theoretical and Observed Frequency Response 
of Bilinear Low Pass Filter; = 0.0 
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FIGURE 2.15c Theoretical and Observed Frequency Response 
of Bilinear Low Pass Filter; = 0.0 
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FIGURE 2.16a 



Theoretical and Observed Frequency Response 
of Bilinear Low Pass Filter; = 0.3 
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FIGURE 2.16b Theoretical and Observed Frequency Response 
of Bilinear Low Pass Filter; b^ = 0.3 
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FIGURE 2.16c 



Theoretical and Observed Frequency Response 
of Bilinear Low Pass Filter; = 0.3 
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FIGURE 2,17a Theoretical and Observed Frequency Response 
of Bilinear Low Pass Filter; = 0.7 
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FIGURE 2.17b 



Theoretical and Observed Frequency Response 
of Bilinear Low Pass Filter; - 0.7 
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2.17c Theoretical and Observed Frequency Response 
of Bilinear Low Pass Filter; = 0.7 
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FIGURE 2.18a Theoretical and Observed Frequency Response 
of Bilinear Z High Pass Filter; bi = -0.3 
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FIGURE 2.18b 



Theoretical and Observed Frequency Response 
of Bilinear Z High Pass Filter; b-j^ = -0.3 
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FIGURE 2.18c 



Theoretical and Observed Frequency Response 
of Bilinear Z High Pass Filter; = -0.3 
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FIGURE 2.19a Theoretical and Observed Frequency Response 
of Bilinear Z High Pass Filter; = -0.7 
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FIGUBE 2.19b Theoretical TIT.T'' 

of Bilinear Z High Pass Ulter, 
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FIGURE 2.19c Theoretical and Observed Frequency Response 
of Bilinear Z High Pass Filter; = -0.7 



68 



TABLE I 



Filter 



Theoretical 




400 kHz f = 1 MHz 
c 



Bilinear Z 



Low Pass 



= 0.0 


0.250 


0.251 


0.246 


0.3 


0.343 


0.338 


0.330 


0.7 


0.444 


0.438 


0.431 


High Pass 








= -0.3 


0.157 


0.153 


0.151 


-0.7 


0.056 


0.061 


0.065 



Standard Z 



Low Pass 



= -0.3 


0.221 


0.239 


0.219 


-0.7 


0.057 


0.072 


0.063 


High Pass 








= 0.3 


0.279 


0.254 


0.2 52 


0.7 


0.443 


0.429 


0.428 
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III. 



SAMPLED ANALOG NON-RECURSIVE FILTER 



A. THEORY 

Since the frequency response of a sampled filter is 
periodic, it may be described by a Fourier series [1], 

H(jo)) = E Dj^ ^ 

k=-oo 



where 



^ ^ do) . 

^ -7t/T 

This expression for D^, is mathematically equivalent to an 
inverse Laplace transform of the transfer function, H(s), 
which is h(t) , evaluted at t = kT, 

Dj^ = h(kT) , 

where h(t) is the impulse response of the filter. By 
using the definition of z, 

00 

H(z) = I h(kT) z^ , 
k=-co 

where h(t) is assumed to be an even function in time, 
therefore it is not a true impulse response. 
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For a delay device with 2N taps available, H(z) can be 
approximated by the series, 

^ k 

H(z) = Z h(kT) 2 ; . 

k=-N 

Figure 3.1 illustrates the implementation of this series. 

The approximation can be made closer to the designed 

response by truncating the series at larger values of N. 

At least 2N storage elements are required. 

For a rectangular bandpass filter response, such as 

that shown in Figure 3.2, repeated with a period of 
2 TT 

, the coefficients of the Fourier series are calcu- 

s T 

lated in Appendix G by the following expression, 

2 

D, = D , = r— sin(2k7TBT) cos (2k — ) 

K *“K KTT Ca)^ 

S 

Since truncating the series is like multiplying the 

infinite series by a window which is unity for -NT t _< NT 

and zero elsewhere, the result of the truncation is to 

convolve the frequency response with the Fourier transform 

of the window function. A rectangular window transforms to 
s in X 

a — - — function which adds sidelobes to the response as 
a result of the convolution. By choosing a different window 
fxinction, such as a Hamming window, these sidelobes can be 
minimized at the expense of greater passband width [8] . To 
use a Hamming window, each coefficient, Dj^, is multiplied 
by a weighting factor, Wj^, where. 



71 




c 

•H 



FIGUl^ 3.1 Implementation of a Non Recursive Filter 
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FIGURE 3.2 Rectangular Bandpass Filter Response 

I 

Wy, = 0.54 + 0.46 cos (it jy) 

Since the input to the non-recursive filter is sampled 
with a finite width pulse of T^, the clock period, the 
spectrum of the input is enclosed within an envelope, E(f) , 

sin(TTfT^) 



This envelope has its first null at f f^/ the clock 
frequency. Though this envelope is imposed by sampling 
rather than the filter algorithm, it can be included within 
the filter response function as a multiplicative factor. 
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B. EXPERIMENT 



1. Configuration 

Reticon Corporation markets a 12-term Tapped Analog 
Device (TAD-12) which provides 12 output taps with the 
following time delay configuration. 




FIGURE 3.3 TAD-12 Delay Configuration 

Since the first tap is only T^ delayed from the input and 
the remainder of the taps are spaced with 2 t^, the undelayed 
input can not be used for the A^^ term, therefore, only 11 
terms of a given tapped delay device can be utilized. 

The filter is implemented by weighting the tap 
outputs proportional to the Fourier series coefficients. 
Analysis performed in Appendix G and the computer program 
described in Appendix H provide design values for the 
tapping resistances, Rj^. Figure 3.4 describes the filter 
set-up. The lower ends of Rj^ are connected either to the 
inverting or non-inverting input of the summer depending 
upon the relative coefficient polarity. Since Dj^ = D_j, , 
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3.4 Non-Recursive Filter Set-Up 
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Rj^ = In order to not load down the taps, a minimum 

value of 10,000 ohms is set for the tapping resistances. 
Appendix H indicates the following expression for Rj^, 

D 

% = (lOK + r) (-5S2) . r , 

l°kl 

where r is the tap output impedance. For the TAD-12, r is 

cibout 5K. Since the resistances depend on the relative 

values of the series coefficients rather than the absolute 

values, gain scaling is accomplished by adjustment of the 

feedback resistance, R^ . For consistent scaling for 

^1 

opposite polarity coefficients, R_ must be equal to R^. 

Due to the TAD- 12 circuitry, the tap outputs are at 

a +5 volt bias level, so the lower end of R^ returns to 

^2 

+5 volts rather than ground. The output of the summer, 
therefore, rides on a 5 volt bias level. 

2 . Observations 

Two sets of coefficients were designed for a Q of 
15. The first set are for fQ = and the other for 

fQ = %fg . Table II lists the coefficients and corresponding 
resistances. These sets of coefficients were measured at 
clock frequencies of 50 and 100 kHz. 

Figures 3.5 through 3.12 illustrate the experimental 
frequency response compared to the theoretical calculations. 
For these data, the measured center frequency, fg, passband 
width, BW, and Q are listed in Table III. 
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TABLE II 



k| 




""Ikl 




""Ikl 


0 


0.03333 


10 K 


0.03333 


10 K 


1 


0.0 


open 


-0.03084 


11.2 


2 


-0.02770 


17 ,0 


0.02419 


15.7 


3 


0.0 


open 


-0.01556 


27.1 


4 


0.00555 


85.0 


0.00773 


59.7 


5 


0.0 


open 


-0.00314 


154.3 







TABLE III 






f 


FFS 


f 


BW 


Q 


c 




o 




50 kHz 


0.25 


6.25 kHz 


3.213 kHz 


1.95 




0.5 


12.5 


3.033 


4.12 


100 kHz 


0.25 


12.5 


6.402 


1.95 




0.5 


25.0 


6.154 


4.06 
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Ma^iiude (noroallzed) 




FIGURE 3.5 Theoretical and Observed Frequency Response; 

“ '•‘^s - ” 50 kHz 
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FIGURE 3.6 Theoretical and Observed Frequency Response; 

= 50 kHz, Magnitude in dB 
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Ma^ltudo (nomallzed) 




FIGURE 3.7 Theoretical and Observed Freqiaency Response; 

= 100 kHz 
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FIGURE 3.8 Theoretical and Observed Frequency Response; 

= -ifg f = 100 kHz, Magnitude in dB 
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FIGURE 3.9 Theoretical and Observed Frequency Response; 

fo = ' fc " 
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FIGURE 3.10 Theoretical and Observed Frequency Response; 

f_ = hfc f = 50 kHz, Magnitude in dB 

O o C 
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FIGURE 3.11 



Theoretical and Observed Frequency Response; 
^o " ^*^s ^ ^c " 
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FIGURE 3.12 Theoretical and Observed Frequency Response; 

- 100 kHz, Magnitude in dB 
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It is noteworthy that Qdid not change significantly when 
the clock frequency was doubled and that the bandwidth 
changed by a small amoimt when the center frequency was 
doxabled by coefficient changes. 

The measured Q was significantly below the designed 
value of 15. Some insight into the reason for this discrep- 
ancy may be found by looking at the impulse response of a 
rectangular bandpass filter. Without rejecting the negative 
time portion of the impulse response, since the non-recursive 
filter uses a 'two-dised' impulse response, one has from the 
inverse Fourier transform of a rectangular window centered 
at plus and minus 



h(t) 



2 "o , 

— COS(0)„t) 



sin (tt 



Wot 

Q 



ir 




1 



This is a continuous function for a non-repeating frequency 
response. For a frequency response periodic with period, 
w_, the impulse response is discrete with the above function 
acting as an envelope. Figure 3.13 illustrates the response. 
In actuality, the sampled nature of the impulse response 
causes the repetitious frequency response, but the duality 
of the Fourier transform ignores the true cause and effect 
relationship. The spacing of samples is. 
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FIGURE 3.13 Two-Sided Impulse Response 

^o _i 

For — = a , there will be a samples per cosine period. 
There are Q divided by pi cosine periods under the main lobe 
of the envelope. Therefore, the n\amber of samples within 
the main lobe is, 



If one uses the criterion that the main lobe must be 
represented completely in order to approximate the impulse 
response, i.e. M samples are required, then some estimate 
of the number of terms required to give a selected Q at a 
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is available from the above 



center frequency of aw^ 
equation. Conversely, if the number of terms is fixed, the 



realizable Q may be estimated given the ratio of tj to w . 

o s 



In this study, M = 11 
lists the expected and 


and a = 0.25 
measured Q's. 

TABLE IV 


and 0.5. 


(JL) 


o 

03 

s 


^expected 


Q 

measured 


0.25 


8.64 


1.95 


0.5 


17.28 


4.1 



The observed Q's did occur in about the expected ratio of 
two to one, but were significantly less than the expected 
values. Analysis in Appendix J indicates that this lowering 
of Q is mostly due to the Hamming weighting used. The 
expected values for this weighting are 2.17 for the center 
frequency to repetition period ratio of 0.25 and 4.35 for 
the 0.5 ratio. These match more closely to the observed 
values. The remaining difference is probably due to the 
need for more impulse response terms than those in the main 
lobe. 

A significant amount of clock noise, up to 40 milli- 
volts, was observed at the output. A low pass filter would 
be desirable at the output to decrease this noise content. 
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Perhaps this low pass filter could be built in conjunction 
with a level shifter to remove the output bias. It was 
further observed that using lead lengths of a few inches 
to connect the coefficient potentiometers to the TAD added 
significantly to the clock pick-up, which severely limits 
the small signal handling capability of the device. Input 
signals beyond 8 volts peak to peak began causing distortion 
in the tap outputs. Therefore, without using a low pass 
filter to suppress the clock noise, the dynamic range of 
this device is less than 54 dB. 
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IV 



SAMPLED ANALOG FOURIER TRANSFORMS 



A. THEORY 

As mentioned in Chapter I, a possible application of 
sampled analog signal processing is to use the combination 
of a recursive comb filter and a non-recursive bandpass 
filter as a sampled analog Fourier transformer. The teeth 
of the comb filter breaks a signal into its Fourier compo- 
nents, while the bandpass filter selects one of the 
components. Figure 4.1 illustrates such a system. 



Comb Filter Bandpass Filter 








Input Spectrum 



Sampled Spectrom 




FIGURE 4.1 Sampled Analog Fourier Transformer 
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B. EXPERIMENT 



For the purpose of demonstration, a bilinear Z- transform 
low pass filter with a narrow tooth width and a non-recursive 
bandpass filter are implemented and cascaded. For this 
example, the comb filter coefficients are chosen as, 

1 , 

1 . 

- 0.7 . 

The design of a generalized analog Fourier transformer 
can be made based on the number of frequency components 
desired, i.e. the number of comb teeth selected. The number 
of desired components, K, determines the frequency of the 
highest order tooth to be selected. 




f 

max 



K f 



s 



Since fo “ ^ / where T = LT for L equal to the number 
Si c 

of delays in the recursive filter delay device and T^ is 
the recursive filter clock period , 



f 

max 




where f 

c 



1 

T * 



Since f^ is the sampling rate, the Nyquist criteria states 
that. 
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I 






so 



f 

max 



< 



% f 



c 



L > 2K 



For a bandpass filter designed for a center frequency to 
repetition period ratio, 







the next highest passband is located at (l-a)w . Since 
the first K teeth should be distinct, the first repeated 
passband should occur at least at a frequency. 



(1 - a) o)g 



^ (K + 1)2tt 
- T 



f 



while. 




These restrictions on a can be summarized. 



a 



< 



1 

K + 2 * 



If the bandpass is tuned by changing its clock frequency, 
its Q will remain constant, so its bandwidth will increase 
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as it is tuned to higher ordered teeth. At the maximum 
value, i.e. ^ fg » the bandwidth will be, 



BW 



max 



K f. 



e f 

s 



The term, B, is a design choice that should be some fraction 
less than vinity. The required Q is, therefore. 



From the equation relating Q, a, and M, the number of taps, 
a guess at the required device length is made. 



M > 




K(K + 2) 
Ba 



Thus, for a possible design to find 100 frequency components, 
and a maximum bandwidth of one tooth repetition period, 
i.e. B = 1 , 



M > IQQ (100 + 2) 

— IT 



3247 . 



This means that the Fourier series is from the -1623^ term 
to the 1623^^ term. 

The non-recursive filter coefficients are selected for 
a center passband frequency of ^ • Figure 4.2 shows the 
resulting transfer function for the cascaded filters plotted 
against the comb response. The comb filter clock was run at 
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FIGURE 4.2 



Observed Frequency Response 
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at 400 kHz, which placed comb teeth at approximately 
4.1 kHz increments. The bandpass filter clock was adjusted 
so that the passband selected the third tooth. The frequency 
of the third tooth maximum was m.easured to be 12.452 kHz. 

The bandpass filter clock was measured to be 49.555 kHz. 

Several aspects of the cascaded response are noteworthy. 
The shape of the passband is determined by the comb filter 
tooth shape, which should remain constant from one tooth to 
the next, whereas the non-recursive bandpass filter bandwidth 
as shown in Chapter III, changes as it is tuned with its 
clock. This use of the comb filter to maintain a constant 
spectrum 'sample' window allows one to keep higher order 
harmonics separated. The other important aspect of the 
system which is noteworthy is the isolation achieved between 
the selected tooth and the unselected teeth. The following 
table lists the first several teeth and the isolation of each 



Tooth 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 



TABLE V 

Isolation 

-23.1 dB 
-27.1 
0.0 
-26.7 
-34.2 
-44.0 
-37.1 
-34.9 
-13. 4 
-44.0 



The low isolation figure for the ninth tooth resulted from 

3 

the repeated response at j f^ . 
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V. CONCLUSIONS 



The Z-transform design methods developed for digital 
filter theory are useful for sampled analog recursive filter 
design. High and lov? pass cut-off frequencies v;ere predict- 
ably determtined by the filter feedback coefficient values. 

These values were adjusted by means of a generalized cali- 
bration curve or by measurement of the unit step response 
of the filter. Som^e deviations from theoretical frequency 
responses were observed. These deviations are apparently 
caused by instrum.entation limitations. 

Non-recursive filters were implemented by truncation of 
a Fourier series expression for the desired repetitive 
frequency response. The center frequency of the bandpass 
filter response passband was tuned at nearly constant Q by 
varying the clock frequency. The pass band was also shifted 
at nearly constant bandv;idth by changes made to the coeffi- 
cients of the filter. The sin X over X envelope caused by 
the finite sampling pulse-width multiplies the desired 
frequency response. This result decreases the signals passing 
through the repeated passbands. 

Several particular frequency components of an input signal 
spectrvim were examined in real time by cascading a recursive 
'comb' filter and a non-recursive bandpass filter to form a 
sampled analog Fourier transformer. An estimate of the 
required delay device lengths for im.plem.entaticn of a sam.pled 
Fourier transformer by the above m.ethod was found for a general 
number of desired com.pone.nts . 
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APPENDIX A 



Derivation of Filter Transfer Function 



For the system diagram illustrated below, 



in' 




FIGURE A-1 



the transfer function, H(z), is defined as. 



H(z) 






From the system diagram it follows that, 

X(z) = E. (z) + (-b,) (X(z) z~h + (-b,)(X(z) z“^) , 

Xli X ^ 
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SO/ 



Ein(z) = X(z) + b3^X(z)z"^ + b2X(z)z"^ . 

The output of the system is, 

Eout^^) = a^XCz) + A^X(z)z”^ + a2X(z)z”^ , 
therefore. 



H(z) = 



Eout^^) X(z) (a^ + a^z”^ + a2Z 

X(z) (1 + h^z~^ + b2z“^) 



u. -1 ^ -2 

aQ + a^z + a2a 
1 + b^z ^ + b2Z ^ 



98 





I 



APPENDIX B 



Recursive Filter Design 

The transfer function of a given second-order sampled 
recursive filter is given by the following expression, 



H(z) = 



-1 -2 
a^ + a^^z + 5 ^ 2 "^ 

1 + b^z~^ + b2z“^ 



By using the definition of z. 



s T "1 (joT 

z = e = e-' , where T is the time delay. 



the frequency response, H(jaj), can be determined: 



H{jo)) 



^0 -I- 

1 + bj^e + b2S 



j2o)T 

J2^ 



/ 



ao+ai(cos wT - j sin wT) + a 2 (cos 2coT - j sin 2(jjT) 
1 + b^ (cos cjT - j sin wT) + b 2 ( cos 2coT - 3 sin 2coT) 



Since the only variation of H(jw) with frequency is due to 

the sine and cosine functions, and since they are periodic 

in w for a fixed T, H(jw) is periodic in frequency. In fact, 

2 TT 2 7T 

the period of the response is (k) ^ ^ ^ (k+ 1 ) , 
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for k = . . . -2 , -1 , 0 , 1 , 2 , . . . . Since the sine and cosine 
teirms repeat themselves exactly, the entire frequency response 

of the filter can be described by looking at a single period, 

2 TT 2 7T • 

say 0 ^ • The frequency, , is called the 

repetition period, w , which is related to the time delay, T. 

o 

As illustrated in Figure 2.6, the frequency response from 
0 < 0 ) < is mirrored for < cu < o) . This periodic 

nature of the filter frequency response, which arises from 
the use of delays, can be used to implement comb filters as 
described above. 

There are several methods for designing the coefficients 
of the transfer function, H(z). Two methods are frequently 
used if the Laplace transfer function, H(s), or the impulse 
response, h(t) , of the desired continuous response is known. 
One is the standard Z-transform procedure. The other is the 
bilinear -Z-transform method. 

For the standard Z-transform approach, h(t) is determined 
for the desired response and then sampled at intervals, T, 
the time delay. 



^sampled 



E h(kT) 6(t-kT) 
k=0 



Its Laplace tranform is. 



^^^sampled^^^ ^ 



E h(kT) e 
k=0 



lOU 



From the definition of z = e 



ST 



00 



H(z) = Z h(kT) z 
k=0 



Z-trcinsform tables are available for those impulse functions 
of interest where a closed form of the series is known. 

The second technique, the bilinear Z-transform, involves 
the use of an approximation for s in terms of z. The 
transfer function, H(s) , is converted directly to H(z) by 
the sxibstitution of. 



The use of this approximation for s leads to a non-linear 
relationship between the frequencies of s = jw and 



2 z - 1 
T z + 1 



z 



e 



jw' T 



_ . _ 2 z - 1 

® T z + 1 



so 




therefore 



= tan isw'T 



or. 
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» 



39 



= In 

T 



Since ~ ~ ' 



-1 

03 ' = ~ tan ^ihoiT) 

TT 



This non-linearity is defined as frequency warping. 

A. STANDARD Z-TRANSFORM DESIGN PROCEDURE 
1. First-Order 

a. Low Pass Filter 

A first order RC low pass filter has the impulse 

response/ 

-0) t 2^ 

h(t) = u)^ e , where ~ 



The standard Z-transforra is, 



<» - 0 ) nT 

H(z) = w E e ^ z"^ 

^ n=0 



w. 



-w T--1 
c z , 

a / — I 

-e. 



1 - 



For a first order low pass filter, 
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Since the value of 3 q is a constant factor which sets the 
filter gain, normalizing the response is equivalent to making 
a normalized transfer fiinction with Sq = 1. Since the 
coefficient is the only term which acts on the frequency 
dependent z, the cut-off frequency, is determined solely 

by b^. From the definition of cut-off frequency. 






= J |H(j03) 1 



2 

max 



1 

1 + b^ (cos w^T - j sin 



2 



1 + 2b, cos CO T + b,'‘ 
1 cl 



The maximum value of H(jco) occurs when 2b^^ cos is 

most negative, i.e. that factor is equal to -2|b^| . 

Therefore , 



1 + 2h^ 



1 

cos CJ^T + 




1 




+ b 



2 

1 



or 



1 + 2b, cos coT + b-^ = 2- 4 lb, I + 2b, ^ 
1 c 1 ' 1 ' 1 



cos co^T = 



1 - 4|b^| + bj^^ 
2hl 
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-1 1 - 4|bj^| + bj^ 

* zf: > 



(JL)^ 






2tt 



cos 



-1 1 - ^l^ll 



+ b. 



2b, 



TABLE VI 



0.1716 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 



0.5000 

0.3524 

0.2207 

0.1573 

0.1150 

0.0831 

0.0574 

0.0357 

0.0168 

0.0000 



For the first order low pass filter, < 0 , 
maximum values will occur whenever 2bj^ cos loT 
when cos toT = 1 , or 



(0 = k(^) , k = ...-2, -1,0, 1,2,. 

= kio^ . 



therefore 
= 2bj^ , or 

• / 
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Likewise, the minimum value of the response occurs when 
ujT = -1 , or , 



0) = (2k+l) (^) 

= (k+%)Wg . 



The ratio of the maximum response to minimum response yields 
an interesting result. 



y 



2 



|H(jw) 

|H(jcj) 



2 

min 



max 



1 

1 + 2|b,| + 

1 

1 - 2|b^| + 

g - |bii)^ 

(1 + IbiD^ 



y = 




Solving this equation for yields. 



b 



1 



1 - y 
1 + y 
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TABLE VII 



U |b^| 



0.0 

0.01 

0.025 

0.05 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 



1.0000 

0.9802 

0.9512 

0.9048 

0.8182 

0.6667 

0.5385 

0.4286 

0.3333 

0.2500 

0.1765 

0.1111 

0.0526 

0.0000 



This result can be of use in experimental determination 
of the effective value of b^. 

b. High Pass Filter 

The low pass filter just described above can 
be transformed into a high pass filter by shifting the 
response by so that the maximvun values occur at 

o 

(k + h) w_ and the minimum values at kw . Changing the 
s s 

polarity of b^^ would accomplish this result since the 



maximum response occurs 


when 


2h^ 


cos 


wT = 


-2|b^| , which 


would make cos wT = -1 


and 


w = 


(k + 


h) Wg 


for b^ > 0 . 
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■i 





if 




The minimxim response would occur when cos uT = 1 , or 
0 ) = . The cut-off frequency, co , is found by the same 

S w 

procedure as for the low pass case: 






||H(jo)) 



2 

max 



1 



- 30 ) T 

1 +b^e 



2 



1 + 2b, cos ( 1 ) T + b' 
1 cl 



From the previous case, since the expressions for the 
cut-off magnitudes are identical. 




^ cos ^ ( 
2tt 



1 - 



4bi + 



2b. 



TABLE VIII 



0.1716 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 



0 ) 

c 

0) 

s 

0.0000 

0.1476 

0.2793 

0.3427 

0.3850 

0.4169 

0.4426 

0.4643 

0.4832 

0.5000 
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For this form of the first order high pass filter, the other 
coefficients are the same as for the low pass case, i.e. 
aQ = 1, aj^ = 0 for a normalized transfer function. 

Another type of high pass filter may be imple- 
mented by the standard Z-transform technique. For a first 
order RC high pass filter the impulse response is, 

h(t) =^^[H(s)] 

"^c*^ 1 

= 6(t) - co^e , where ^ • 



The Z-transform of this impulse response is. 



00 —CO nT _ 

H(z) = Z (6(nT) - WcS ^ ) z ^ 

n=0 



, -1 
1 - “c - e 2 

, -“c^ -1 

1 - e z 



The transfer coefficients are. 




The corner frequency is a function of both a^^ and b^^. 
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2. Second-Order 



a. Low Pass Filter 

For a second-order LC low pass filter, the 
impulse response is 



0) ^ 

h(t) [ — - 9 .-. ■] 

c 



w sin w t 
c c 



where w 

c 




The resulting Z- transform is, 



H(z) = E w sin w nT z 

O c c 



oj sin 0 ) T z 
_ c c 

“ -1 -2 ' 

1-2 cos (jj T z + z 
c 

For a normalized transfer function, Uq = a 2 = 0 , bj|^ = 1 , 
and, 

= 1 f 

b, = -2 cos w T . 
i c 

From the expression above, the relationship between b^^^ and 
the corner frequency, oj , can be stated. 




1 -1 
^ cos i-hh^) 
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TABLE IX 



-2.0 


0.0000 


-1.8 


0.0718 


-1.6 


0.1024 


-1.4 


0.1266 


-1.2 


0.1476 


-1.0 


0.1667 


-0.8 


0.1845 


-0.6 


0.2015 


-0.4 


0.2180 


-0.2 


0.2341 


0.0 


0.2500 


0.2 


0.2659 


0.4 


0.2820 


0.6 


0.2985 


0.8 


0.3155 


1.0 


0.3333 


1.2 


0.3524 


1.4 


0.4734 


1.6 


0.3976 


1.8 


0.4282 


2.0 


0.5000 



b. High Pass Filter 

The impulse response for a second order LC high- 
pass filter is. 



h(t) = [-~ 2 ^ — 2 ^ = 6(t) - 0 )^ sin u^t . 

s'^ + w 

c 



This transforms into. 



H(z) = 1 (6(nT) - 0 ) sin no) T) z 

n=0 ^ ^ 



-1 -2 

1 - (2 cos 0 ) T + w_sin w T) z + z 

WWW 

_ __ 

1-2 cos (j T z + z 
c 
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1 



where. 



ao = 1 , 



a, = -2 cos 03 T + (jj_ sin o) T , 

I ' 



32 = 1 , 



b, = -2 cos 03 T , 
1 c ' 



b2 = 1 * 



The cut-off frequency for this filter will be dependent on 
both and b^^. 



B. BILINEAR Z-TRANSFORM DESIGN PROCEDURE 
1 . First-Order 

a. Low Pass Filter 

The transfer function of a first-order RC low 
pass filter is, 



03 



H(s) = 



S + 03 



2 z - 1 



By sxabstitution of the bilinear approximation, s = ^ ^ / 
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H(Z) 



“c 

2 z - 1 ^ 

T z + 1 “c 



(1) T(z + 1 ) 



2z “ 2 + 03 (z + 1) 
c 



co^T + o)^T z 



-1 



2 + (0^ + (o)^ “ 2) z 



-1 



^-1 



2+cOg • 2+0)^ 

“2 

^ 0) T + 2 ^ 



cOcT 



Normalizing the transfer function by 75— r yields, 



H(z) = 



1 z-1 

“o'^ ■ 2 -1 



where, 



Sq ^ ' 



a^ = 1 , 



h - 

1 (o_T + 2 



The bilinear Z-transform warps the frequency response, so 
the actual cut-off frequency, o)_,, is , 

W 



“c' 



I tan"^(%co^T) 
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From the expression for in terms of 



-2 

"c T *bj^ - 1* ' 



SOi 






-2 ^ - 1/^1 ''' 

_ tan (^-rr) . 



or. 



03 , 

C ' 

03 



-1 



TABLE X 



b. 


* 


1 


“s 


-1.0 


0.0000 


-0.9 


0.0167 


-0.8 


0.0352 


-0.7 


0.0556 


-0.6 


0.0780 


-0.5 


0.1024 


-0.4 


0.1289 


-0.3 


0.1572 


-0.2 


0.1872 


-0.1 


0.2183 


0.0 


0.2500 


0.1 


0.2817 


0.2 


0.3128 


0.3 


0.3428 


0.4 


0.3711 


0.5 


0.3976 


0.6 


0.4220 


0.7 


0.4444 


0.8 


0.4648 


0.9 


0.4833 


1.0 


0.5000 
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b. High Pass Filter 



For a first-order RC high pass filter the 
transfer function is, 



H(s) 



s 

S + Cl) ' 
c 



which transforms to. 



H(z) 



2 z - 1 
T z + 1 
2 z - 1 
T z + 1 



+ 



0) 



c 



z - 1 

- 1 + (z + 1) 



-1 



co^T + 2 



0) T + 2 
c 



Wc"^ 2 

1 + — — — ^ Z 



“o’f + 2 



For the normalized transfer function. 



a^ = -1 , 

0) T - 2 
c 

^1 ^ w T + 2 * 

c 
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It is noteworthy that the only change from a high pass filter 

to a low pass filter, other than the gain normalization 

factor, is the polarity of a^^. Since the expression for 

in terms of the design cut-off frequency, w , is the same 

as for the low pass filter, the same relationship between 

b, and the actual cut-off frequency, lo , , as found for the 
J- c 

low pass case exists. 

2 . Second-Order 

a. Low Pass Filter 

The transfer function for a second-order LC low 
pass filter is. 



H(s) = 



2 — — 2 • where = (LC) ^ 

s'^ + 0) 



The bilinear Z-transform for this function is. 



H(z) = 



w. 



f2 z - 1 ,2 . 2 

T z + 1 ^'c 



w^^T^(z^+z+l) 



4z^ - 8z + 4 + CO ^T^(z^+2z+l) 

o 



2 m 2 o 2_2 2^2 

03c T . T co^ T ^_2 

2„2 , , 2„2 . . ^ O 



W(, 'T +4 T‘ + 4 



-T +4 



- 8 . .2 
1 + — z + z 

CO ^T^ + 4 
c 
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where, for a normalized transfer function, 



aQ - 1 / 
— 2 / 
a2 = 1 , 



b 



1 




- 8 



+ 4 



b 



2 



1 . 



By using the expression for b^^ to solve for ^ in terms of 
bj^, one finds. 



0) 



c 



2 2 + bj^ 

T '2 - b ‘ 



Because of the frequency warping. 



2 + b. 1 



2 4 - “ 1 / 

“c- = T <2^ 



1)- 



or 



0) , 

c ' 

03 



1 -1,^ 



It is interesting that this design table is the same as 
Table IX, the standard Z-transform second-order low pass 
filter. 
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TABLE XI 





o)c. 


b. 


1 


‘^s 



- 2.0 


0.0000 


- 1.8 


0.0718 


- 1.6 


0.1024 


- 1.4 


0.1266 


- 1.2 


0.1476 


- 1.0 


0.1667 


- 0.8 


0.1845 


- 0.6 


0.2015 


- 0.4 


0.2180 


- 0.2 


0.2341 


0.0 


0.2500 


0.2 


0.2659 


0.4 


0.2820 


0.6 


0.2985 


0.8 


0.3155 


1.0 


0.3333 


1.2 


0.3524 


1.4 


0.3734 


1.6 


0.3976 


1.8 


0.4282 


2.0 


0.5000 



b. High Pass Filter 

A second-order LC high pass filter has the 
following transfer function, 




9 



which is transformed to, 
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H(z) 



f 2 z - 1, 2 
'■T z + 1^ 



,2 z - 1, 2 , 2 

TTT^ + 



4 (z - 2z + 1) 



4z^ - 8z + 4 4- (z^+2z+l) 



3z 



-1 



4z 



-2 



2 2 

+ 4 



2 2 2 2 

4- 4 + 4 

c c 



^ ® -1 -2 
1 4- z 4- Z 



+ 4 



For a normalized transfer function, 



aQ - 1 / 



a^^ — “2 / 



So = 1 / 



b, = 



2 2 

2w T - 8 
c 

2 2 

0) + 4 

c 



= 1 . 



It is noteworthy that this transfer function has the same 
coefficients as the low pass filter with the exception of 
the sign of a^^. The corner frequency is dependent on b^^ 
in the same manner as for the low pass case. It is now 
apparent that the polarity of a^^ is the means for converting 
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high pass bilinear Z-transforras to low pass and vice versa. 
In addition, the aQ, a^^, and a2 values are always integral 
multiples of one another. 

C. SUMMARY 

The following table summarizes the design procedures 
for the standard Z-transform and the bilinear Z-transform 
for first order filters. 
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TABLE XII 



Design Method 



Standard Z 



Bilinear Z 



Filter Type Coefficients 

a^ 



Response 

Period 



Lew Pass 



— < J 
03s 4 



r\ -03 T 

c -e 



IcM Pass 



cOs 4 



-wT -to T 
, c c 

1 “ to -e e 

c 



High Pass 

— < i 

tOg 4 



-to T -to T 

1 -to -e ^ -e 
c 



Hi^ Pass 

to 4 
s 



to 



-to T 
c 






Low Pass 




Low Pass 




Hi^ Pass 




Hi^ Pass 




to T 
c 




to T 
c 


to T-2 
c 




to T + 
c 


2 


to T + 2 
c 


to T + 2 
c 


to T 
c 




CO T 

c 


to T- 2 
c 




to T + 
c 


2 


to T + 2 
c 


to T + 2 
c 


YY 


2 




-2 


<N 

1 

3 


nn 


to T + 
c 


2 


to T + 2 
c 




2 




-2 


<a^T-2 




to T + 
c 


2 


to T + 2 
c 


to T+ '2 
c 


AA 
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APPENDIX C 



Derivation of Calibration Curve 
For Recursive Filter Coefficients 

The coefficients of the recursive filter transfer 
function are produced by potentiometers. These potentiometers 
are loaded by the input impedance of the summing circuit and 
are driven by the output of the SAD-100 driving board as 
illustrated below. 




FIGURE C-1 Coefficient Circuit 

It is assumed that the summer op amp is operated as a unity 
gain, inverting device. Therefore, the voltage on the tap 
of the pot is The voltage on the top of the poten- 

tiometer is -av. . The coefficient, b, is defined as the 

in 

ratio of to From the voltage divider equation. 
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-^out = -^^in + 



R| 

R I IR^ ' 



therefore. 



-V 



b = 



out 



aR_R 

cl 

R^ + R 

ci 



^in 






R^ + R 

O. 



= a 



''a'' 



R^(Rp-R^) ^ RpR 



The limit of the coefficient, b, as R, the loading resistance, 
gets arbitrarily large is. 



lim b 

R ->oo 



lim a 
R -»<>= 



R^R 

cL 



R (R - R ) + R R 
a p a p 




This is the expected value for a non-loaded voltage divider. 

Since R is finite in this study, its effect must be 

considered. For the actual hardware used, R = lOK, 

P 

R = 30. 9K, and a = 1. The following calibration curve was 
calculated. 
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0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 



TABLE XIII 



b 


^observed 


0.000 


0.000 


0.097 


0.095 


0.190 


0.180 


0.281 


0.286 


0.371 


0.369 


0.463 


0.454 


0.557 


0.546 


0.655 


0.639 


0.761 


0.734 


0.875 


0.858 


1.000 


1.000 
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APPENDIX D 



DETERMINATION OF TRANSFER FUNCTION 
COEFFICIENTS FROM UNIT STEP RESPONSE 

A particular filter configuration can be characterized 
by a transfer function, H(z), where 

-1 -2 
a^ + a, z + a« z 

H(z) = — , 

l + bj^z-^ + b2Z^ 

s T 

for z defined as e and T a time delay. Since a transfer 
fianction relates the output of a system to its input. 



H(z) 



E . (z) 
out 

W. vT5 

in 



i 



or. 






^ -1 ^ -2 
a. + a, z + a^z 

= E. (z)- 

1 + b,z ^ + b.,z ^ 



Multiplying the above equation by the denominator of H(z) 
yields , 



E . (z) + b, E ^(z) z*"^ + bo E . (z) z”^ 

out 1 out 2 out 






-*? 
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so, 



-1 -2 -1 
E . (z) = a»E . + a,E. z + a_E. 2 - b,E .z 

out 0 in 1 in 2 in 1 out 



- b^E .2 
2 out 



-2 



Taking the inverse Z-transform yields the following time 
expression, 



e . (n) = a^e. (n) 
out 0 m ' ' 



f a,e. ^n-i; 
1 in 



2 in 






1 out 






where. 



e . (i) = 
out 



e . ( iT ) 
out ' 



Assuming zero inputs and outputs for n £ 0 and a unit step 
input, i.e. e^^(n) = 1 for all positive n, the output of 
the filter will be the following sequence for the first 
five values of n. 



Table XIV 



n 



e . (n) 

in 



e . (n) 
out 



0 

1 

2 

3 

4 



0 

1 

1 

1 

1 



% * H - ’= 1^0 

■ ‘>l'V®r’=l®0>‘ ‘>2^0 

“^ 2^0 ^ ”*^2 

-b^ (a^+ai-b^a^) -b2ao-b2 (a^+aj^-bj^a^) ) 



-b2(ao+a^+a2-bj^(aj+aj^-bj^aj)-b2a„) . 



Defining the difference between successive steps in the 
output wavefore as “ ®out^^”^^ ' following 

relations are found. 



di = ao - 0 = ao , 



^2 = ^ 0 ’"^ r ^ l ^ 0"^0 = ^1 " ^ 1^0 ' 






= a2 - b^(a^-b^aQ) - b2aQ , 
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similarly, 



The relationship for higher differences is, 



^ ~^1 \-l ~ ^2 \-2 ' k = 4,5,6,... . 



It is clear that the term, subtracted out in 

each difference calculation after k = 4, thus the higher 
differences are a function of only b 2 » and previous 

differences. Solving this relation for b^^ yields. 



b 



1 




-b 

“d 



2^k-2 

k-1 



for k > 4, 



so , 



^=1 = 



-34 - 



^ 2^2 



From the higher difference relation. 



35 = -b^dj 



Solving for b 2 and substituting the above relation for b^^, 
one finds. 
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^2 



■^- 4 ~ - '^ 3 ^ 
^3' - ‘^2‘^4 



Substituting this expression into the relation for 
yields , 



. _ ^ 2^5 ■ ^ 3^4 

^1 2 

^3 - ^ 2^4 



The a coefficients are now easily determined. 



a 



0 




/ 



= d 2 + 



d 2 + dj. 



*^2S °3'^4 

■^ 3 " - ‘^ 2'^4 



d,d„dc + " d,d-d. 

125 23 134 



^ 2 % 



dj - dj d^ 



aj = d3 + bjdj + bjdj^ 



<^ 3 ^ - <^2 '^4 



All coefficients are now determined solely as functions of 
the differences between successive steps in the output. 
Higher order differences could be used to determine b^ 
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and b2/ but using and d^ minimizes the differences required 
to determine the five coefficients. 

In the case of a first-order transfer function, i.e. 
a2 = b2 “ 0 , a simplified set of equations would result. 



= d. 



a- = 



dj + bj^dj + bjd^ 



0 = d^ + ^ , 



so. 



b 



1 




/ 



a^^ = d2 + 



^2 - ‘> 1*^3 



Furthermore, since d, = 



■‘’iVx ■*’2^2 1=2 = 



>=1 = -d 



k-1 



This is a useful result for experimental determination of 
the transfer function coefficients through measurement of 
the step response since b^ may be determined for several 
values of k and then averaged to reduce error. 



APPENDIX E 



PROGRAMS TO CALCULATE THE RECURSIVE FILTER 
COEFFICIENTS FROM THE UNIT STEP RESPONSE 

The formulae resulting from the derivation in Appendix 
D lend themselves to a programmable algorithm for solution. 
The programmable Tektronix 31 calculator is used to imple- 
ment the algorithm in the laboratory. The first five step 
differences from a unit step input are entered in order. 

The calculator then determines and prints the five second- 
order coefficients of the transfer function, a^, a^^, 
b^, and h 2 - 

The basic instructions for using the program are listed 
below: 

1. Turn on calculator and press "clear". 

2. Insert tape cartridge containing filter 
coefficient programs. 

3. Press "reset" and transfer the desired program by 
pressing "from tape" and the number of the track 
the program is stored on (0 for second-order, 

2 for first-order) . 

4. When "busy" light extinguishes, enter d^^ on 
keyboard and press "cont". (Note: the calculator 
should print out d^^ . ) 

5. Repeat step 4 for ^2 through d^. The calculator 
will finish execution of the program after "cont" 
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is pressed for (d^ for first-order 



program) 



After the coefficients are calculated and printed, the 
program returns to its beginning and is ready to read a 
new d^. 

Two limitations of the second order program are note- 
worthy. First, the determination of all coefficients 
except a^ is dependent on b 2 which is, in turn, dependent 
on the higher order differences. For the coefficients of 
interest in this study the differences became steadily 
smaller, therefore, the higher ordered ones contain larger 
percentage error. This error then propagates through the 
other coefficients. Careful measurement of the step 
response can minimize this problem. 

The other limitation of the second-order program concerns 
2 

the factor d^ “ *^2^4* ^ first-order filter. 



d 



3 



d 



4 



b 



1 



d 



d 



2 



3 



or. 



d 



2 




3 



thus 



d 



2 




3 
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% 



This factor, however, appears in the denominator of the 
expressions for all the coefficients except Even if 

experimental error kept the denominator from equaling zero, 
the answer would be useless. For this reason the first- 
order program was written. As noted above the first-order 
program requires four differences to be entered. The 
coefficient, b^^, is calculated by the following weighted 
average. 



V, = itl 

*’1 - - 4 dj - 4 dj • 

The average weighting is a crude compensation for the 
greater percentage error in the d^ measurement. 
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APPENDIX F 



COMPUTER GENERATION OF THEORETICAL 
RECURSIVE FILTER FREQUENCY RESPONSE 



For a given recursive filter transfer function, H(z), 
where, 



H(z) 



^ -2 

ao + ai z + a^ z 

ZT 

1 + z + b 2 z 



/ 



the frequency response, HCjw), may be calculated by using 
the definition. 



z = e 



sT 



j (joT 

= e-^ = cos 



(coT) + j sin (o)T) 



where T is the time delay. By substitution. 



a» + a, (cos coT - j sin wT) + a„(cos 2wT - j sin 2coT) 

TT / -I \ _ V -L ^ 

“ 1 + bj^Ccos u)T - j sin coT) + b 2 (cos 2uT - j sin 2uT) 

a^ + a^^cos ojT + a2COS 2 ojT - jCa^^sin coT + a2sin 2u)T) 

1 + bj^cos o)T + b 2 COS 2coT - j (b^^sin coT + b 2 sin 2ooT) 

I 

For a given frequency, the magnitude and phase of the filter 
response can be determined. 



2 

H(jw)|= [ReH(ju))] + [ImH(joj)] 



/ II / • \ m “1 rImH ( joj) I 

^H(:co)=Tan 
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A computer program was written to calculate, given the 
transfer function coefficients, the response m.agnitude and 
phase over a specified number of steps of a specified 
interval of frequency. The results were then plotted on 
the Calcomp plotter. Experimental magnitude versus fre- 
quency measurements were read in, frequency normalized, and 
super-imposed over the theoretical curves. 

A sample data deck follows: 



1. Line Printer Plot Characters 




* X Y 



1 2 3 4 5 



2. Transfer Function Coefficients 




0 



-0.7 



0 



/ 



1 



10 11 



20 21 



30 31 



40 41 



All nximbers F10.5, in order: a^, a^, a 2 r b^, b 2 



3. Experimental Data Head Card 




13 F10.5 



# Of points Response Repetition Period (kHz) 



4. Experimental Data Points 



0.387 0.001 

/ 

I 1 ... 10 11 ... 20 

' F10.5 F10.5 

Frequency (kHz) Magnitude (normalized) 

The program will print on the line printer the coeffi- 
cients, the data points, a theoretical magnitude versus 
frequency plot, a theoretical phase versus frequence plot, 
and a message stating the status of the Calcomp plot, if 
requested. A labeling statement is used to title the 
Calcomp plot with the coefficients of the transfer function. 
The title declaration cards need to be changed as different 
coefficients are used for the title to coincide with the 
response determining coefficients. 
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APPENDIX G 



Non-Recursive Filter Design 



The design of a non-recursive filter is based upon the 
periodic nature of a sampled filter frequency response. 

The response may be described by a Fourier series, 



H(jw) = E D e 
k=-~ ^ 



where 



tt/T 

i / H(ju) diD . 

^ -1I/T 

The particular filter response of interest in this study is 
a bandpass filter with rectangular passpands. 



-f 



H(f) 

-1 



BP 



j- 



FIGURE G-1 Rectangular Passband Filter 
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which can be described with the parameters: response 

repetition period, passband center frequency, 

and pass bandwidth, BP. The filter quality factor, Q, 
is defined as. 




and the fraction of the response repetition period which 
corresponds to the center frequency is defined as. 




The filter description is complete by specifying f , Q, and 

s 

FFS since the center frequency and bandwidth can be calcu- 
lated from tne above three parameters . The repetition 
period is , 




or Wg = radians/second . 



For a rectangular bandpass filter 



1 



I f - £o I 1 <sBP 



H(jlij) 



0 elsewhere 



therefore 
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D = ^ 
2tt 



2tt (fQ+%BP) 

/ 1 . 
2TT(fQ-}sBP) 



2 -jk(2TrfQT) 

— sin(kTT BP T) e 
kTT 



= ^ sin(kTT BP T) (cos 2-iTfQT - j sin 2TrfQT) . 

Since only real terms may be implemented, the series coeffi- 
cients are, 



Dk = sin(kiT BP T) cos(2TrfQT) . 



The series expansion has infinite terms; however, the non- 
recursive filter has a finite number of taps available to 
approximate the series. 



N 

H(jw) = Z D, e 
k=-N 



jktoT 



N 

Z D, (cos(ktoT) - j sin(kwT)) 
k=-N ^ 



By using the definition of. 



sT j wT 
z = e = e-^ , 
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N 

H(z) = E z 

k=-N ^ 



k 



Since the Z-transform must contain only non-positive powers 
of z, i.e. no anticipated signal values only delayed values. 



This relation states that the transfer function will contain 
a constant delay of NT from input to output in addition to 
the response phase. A possible implementation of this function 
is illustrated in Figure 3.1. From the expression derived 
for the bandpass filter. 



it can be seen that this is an even function of k, so. 



In addition, 

lim = 2-BP-T- (1) • (1) = 2-BP-T , 

k-> 0 

this can be expressed in terms of the design parameters, Q 
and FFS, 



H(z) . z 




k-N 




sin(kTT BP T) 
kir BP T 



cos ( 2k7rf qT) 



D, 



0 
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Window weighting functions are used to modify the 
coefficients in order to improve the filter characteristics. 
A Hamming window is frequently used where the weighted 
coefficient, is, 

Aj, = Dj^ (0.54 + 0.46 COS it)) 

The expression multiplying Dj^ is called the window weighting 
function. 
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APPENDIX H 



Computer Calculation of Coefficients and 
Tapping Resistances for the Non-Recursive 

Filter 



From the analysis performed in Appendix G for a bandpass 
filter. 



\ = D_j^ = 2 cos(kY) , -N < k < N , 



where. 



X = ir , BP = 
^s f = 



pass bandwidth 

response repetition period. 



Y = ^ , fQ 

s 



= passband center frequency. 



Q = ^ 



0 

BP 



FFS = 



For a Hamming weighting window. 



Aj^ = A_j^ = Dj^(0.54 + 0.46 cos (- tt) ) 



For the TAD-12 analog delay device, the coefficients are 
implemented by input resistances to the operational amplifier 
summing circuit, as shown in Figure H-1. 
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FIGURE H-1 Coefficient Implementation 



In order to minimize the loading effect on the tap, the 
tapping resistors, Rj^, are greater than or equal to 10 K ohms. 
The contribution of a single tap, Vj^, to the output, 
is the product of Vj^ and the coefficient, Aj^, or. 




V 



out 



V,. 



r + R, 



This expression indicates that Rj^ is inversely proportional 
to Aj,, therefore the largest coefficient will correspond to 
the minimum allowable resistance, lOK, so, 

A = 

max r + lOK 
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Solving this equation for and substituting that result 
into the Aj, expression yields, 

r + 




or. 




(r + lOK) 



r 



For the TAD-12, the output impedance for a tap, r, is about 
5K. 

A computer program was written to calculate the weighted 

and vinweighted coefficients and their respective tapping 

-1 

resistances for a given T = f , Q, FFS, and a pair of N. 

s 

These parameters are inserted directly into the program as 
follows , 



1. T 




= 0.001 



1 .... 6 7 

( 



13 



2 . FFS 




FFS = 0.25 



1 ... 6 7 ... 14 
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i 



3. 



Q 



Q = 15. 

1 ... 6 7 . . .11 



4. Truncation points, NMl and NM2 

// NMl = 5 

1 ... 6 7 . . . 

/ NM2 = 11 

1 ... 6 7 . . . 

Note: The program calculates results for, 



and. 



H(z) 



NMl 

I 

k=-NMl 






^k-NMl 



/ 



H(z) 



NM2 

Z 

k=-NM2 



k-NM2 

z 



5. Maximum passband magnitude, CA 

/ CA = 1. 

] 

1...67. . . 



Note: CA was set to one in the derivation in Appendix 

G, but the program is more general in this 



respect. 



APPENDIX I 



Computer Calculation of Non-Recursive Filter 
Frequency Response 



Since the transfer function for a non-recursive filter 



IS, 



N ) 

H(z) = E A, 2 
k=-N ^ 



u jwT 

, where z = e-' , 



it is possible to calculate the magnitude and phase of 
H(jo)), where, 



H(jto) = 



N 

E A. e 
k=-N 



jkuT 



N 

E A, (cos(kwT) + j sin(kuT)) 
k=-N ^ 



N 



= Aq + E [ (Aj^+A_j^) cos kojT + j (Aj^-A_j^) sin ktoT] 
k=l 



The real part of H(jto) is clearly the cosine terms, 



N 

Re H(jto) = Aq + E cos(ktjT) , 

k=l 



while the imaginary part of H(jco) is the sine term portion. 



N 

Im H(jo)) = E (A, -A . ) sin(ktoT) 
k=l ^ ^ 
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The magnitude and phase of H(j ) are readily calculated. 



lH(jo))l = [Re H(jw)]^ + [ImH(jw)]^ 






tan 



-1 



, Im H ( jco) 1 
^Re H ( jtjj) ^ 



A program was written to calculate the magnitude and 

2 7T 

phase of H(ju) from 0 = ju = — , plot the result on the 

line printer and Calcomp letter for a given set of Aj,. In 

addition, experimental data points may be read and plotted 

over the theoretical curves. 

The above analysis assumes perfect impulse sampling, 

however, actual devices will have some sampling envelope, 
s in X 

usually a function, which will multiply the magnitude 

X 

of the response. Provision is made in the program to 
introduce this envelope. 

The following sample deck will illustrate the use of the 
program. 

1. Line Printer Plot Characters 



. * X Y 
1 2 3 4 5 
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2. NMAX, the number of coefficients 



NMAX = 2N + 1 



1 1 

12 3 
(13) 



3. Coefficients 



< 






'' -0.0235 

1 ... 10 

(F10.5) 



4 . Data Header Card 

NPTS, the number of data points (multiple of 30) 

FCK, the clock frequency 

XX , the value of sampling envelope at 
magnitude maximum frequency 




50. 0.9845 

... 13 14 ... 23 

FCK XX 

(F10.5) (Flo. 5) 
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5. Data Points 



FX, observed frequency 

HX, normalized magnitude 

6.735 0.033 

1 ... 10 11 ... 20 

FX HX 

I (Flo. 5) (F10.5) 

Note: If dummy cards are required to make NPTS 

a multiple of 30, use FX = 999., HX = 1. 

When the Calcomp plotter is used the "call draw" cards 
are included and the title declaration cards are modified 
to title the plot with the proper number of terms and FFS 
value. The program will always print the number of coeffi- 
cients, the coefficients, and the data points so that proper 
data entry can be checked. 
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APPENDIX J 



Effects of Hamming Weighting 

The use of a weighting factor, Wj^, in the time domain 
is equivalent to convolving the frequency response of the 
desired filter with the Fourier transform of the weighting 
fxanction. The Hamming weighting function, 

Wj,(t) = (0.54 + 0.46 cos(j^ kT) ) , 

which has a Fourier transform, Wj^(jw) , 

0 ! 03 

Wj^(jw) = 0.54 6(jw) + 0.23(6(jw+^) + 6(jw--^)), 



where. 



‘"s = 



2tt 

T 



This result can be illustrated as three impulses as in 
Figure J-1. 



0.54a H(w) 



V’ 



0.25 



-W 



2N 



T 

-J w 



o 

2N 



FIGURE J-1 Fourier Transform of Hamming Function 
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When a bandpass filter is convolved with the impulses, its 
bandwidth is increased by an amount proportional to the 
spacing between the impulses. This added bandwidth will be 
a fraction, y, of the Hamming function bandwidth, 

w -W CO 

added BW = y( 5 | - ^) = -y-j£ . 

The value of gamma is estimated by the following procedure. 
If the values of the impulses were equal, then the -3 dB 
points of the convolved response would be shifted, roughly, 
by spacing of the impulses, ^ , as illustrated below. 




FIGURE J-2 Equal Impulses 



The added bandwidth would be equal to the spacing between 
the impulses, — , or , 



Y = 1 . 
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For the case where the outboard impulses are zero amplitude, 
the response is just convolved with a single impulse and is 
exactly reproduced, so zero bandwidth is added. In this 
case, gamma is zero. A relationship is suggested by this 
bracketing of gamma, 

_ amplitude of outboard impulses 
amplitude of inboard impulse ' 

thus for a Hamming function. 



0.23 

0 . 



0.426 . 



The expected Q for a bandpass filter is. 



Q = 



BW 



Since the Hamming weighting adds some to the bandwidth, B, 
the Q is expected to be less. 



Q' 



CO 



0 



BW + B 



< Q 



By dividing the expression for Q' by that for Q yields 
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Q' 



= Q 
= Q 



= Q 



BW 



BW 


+ B 




BW 




YW_ 


BW 


^ N 




aN 


aN 


+ yQ 



where 



a = 



OJ, 






It is noteworthy that the expression for Q' indicates that 
the limit of Q' as N gets larger is Q. 

If one uses the main lobe criterion for expressing the 
impulse response by specifying the responses within the first 
nulls, one finds that, 



Q = (2N + 1) air 



Substituting this result into the expression for Q' yields, 

n* = ^ (2N^ +N)ir 

^ (2irY + l)N + Y’n’ 

For the case of N = 5 and y = 0.43, 

Q' = a (8.69) 

so for a = 0.25 and 0.5, Q' will be 2.17 and 4.35, respectively. 
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